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Abstract 



Prediction and control of cancer invasion is a vital problem in medical science. This 
paper proposes a modern geometric Ricci-flow and entropy based model for control of 
avascular multicellular tumor spheroid growth and decay. As a tumor growth/decay 
control tool, a monoclonal antibody therapy is proposed. 
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1 Introduction 

Cancer is one of the main causes of morbidity and mortality in the world. There are 
several different stages in the growth of a tumor before it becomes so large that it causes 
the patient to die or reduces permanently their quality of life. Developed countries are 
investing large sums of money into cancer research in order to find cures and improve 
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existing treatments. In comparison to molecular biology, cell biology, and drug delivery 
research, mathematics has so far contributed relatively little to the area pQ. 

On the other hand, consider the vital problem of prediction and control/prevention 
of some natural disaster (e.g., a hurricane). The role of science in dealing with a phe- 
nomenon/treat like this can be depicted as the following OUPC feedback-loop (see [2]): 



with the following four components /phases: 

1. Observation, i.e., monitoring a phenomenon in case, using experimental sensing/measuring 
methods (e.g., orbital satellite imaging). This phase produces measurement data 
that could be fitted as graphs of analytical functions. 

2. Understanding, in the form of geometric pattern recognition, i.e., recognizing the 
turbulent patterns of spatio-temporal chaotic behavior of the approaching hurricane, 
in terms of geometric objects (e.g., tensor- and spinor-fields) . This phase recognizes 
the observation graphs as cross-sections of some jet bundles, thus representing the 
validity criterion for the observation phase. 

3. Prediction: when, where and how will the hurricane strike? 

Now, common, inductive approach here means fitting a statistical model into em- 
pirical satellite data. However, we know that this works only for a very short time 
in the future, as extrapolation is not a valid predictive procedure, even if (adaptive) 
extended Kalman filter is used. Instead, we suggest a deductive approach of fitting 
some data into a well-defined dynamical model. This means formulating a dynami- 
cal system on configuration and phase-space manifolds, which incorporates all pre- 
viously recognized turbulent patterns of the hurricane's spatio-temporal behavior. 
Once a valid dynamical model is formulated, the necessary empirical satellite data 
would include system parameters, initial and boundary conditions. So, this would be 
a pattern-driven modelling of the hurricane, rather than blind data-driven statistical 
modelling. This phase is the validity criterion for the understanding phase. 

4. Control: this is the final stage of manipulating the hurricane to prevent the destruc- 
tion. If we have already formulated a valid geometric-pattern-based dynamical 
model, this task can be relatively easily accomplished, as 




U nder standing 



Prediction 



Control 

i 



Control System 



Dynamical System + Controller. 



So, here the problem is to design a feedback controller/compensator for the dynam- 
ical model. This phase is the validity criterion for the prediction phase. 
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Since there are three distinct stages to cancer development: avascular, vascular, and 
metastatic - researchers often concentrate their efforts on answering specific OUPC- 
related questions on each of these stages (see |3J). In particular, as some tumor cell 
lines grown in vitro form 3-dimensional (3D, for short) spherical aggregates, the relative 
cheapness and ease of in vitro experiments in comparison to animal experiments has made 
3D multicellular tumor spheroids (MTS, see Figure [TJ as well as e.g., Figure 6 in [3]) very 
popular in vitro model system of avascular tumoral . They are used to study how local 
micro-environments affect cellular growth/decay, viability, and therapeutic response (see 
[5]). MTS are often combined with 3D medical imaging |9], 4D confocal imaging |10| . 
3D video holography through living tissual [11] and 3D metabolic imaging imaging bio- 
luminescenc^l [6]. MTS provide, allowing strictly controlled nutritional and mechanical 
conditions, excellent experimental patterns to test the validity of the proposed mathemat- 
ical models of tumor growth/decay [8]. 

A number of mathematical models of avascular tumor growth inside the MTS were 
reviewed in [3]. These were generally divided into continuum cell population models 
described by diffusion partial differential equations (PDEs) of continuum mechanics [2|ll6j 
combined with chemical kinetics, and discrete cell population models described by ordinary 
differential equations (ODEs). Besides, in many cell population models it is possible to 
empirically demonstrate the presence of attractors that operate starting from different 
initial conditions |12} [T3] . 

A general model of multi-phase tumor growth (inside the MTS) is given in [3] by the 

1 In vitro cultivation of tumor cells as multicellular tumor spheroids (MTS) has greatly contributed to 
the understanding of the role of the cellular micro-environment in tumor biology (for review see [SJ 0])- 
These spherical cell aggregates mimic avascular tumor stages or micro-metastases in many aspects and 
have been studied intensively as an experimental model reflecting an in vivo-like micro-milieu with 3D 
metabolic gradients. With increasing size, most MCTS not only exhibit proliferation gradients from the 
periphery towards the center but they also develop a spheroid type-specific nutrient supply pattern, such 
as radial oxygen partial pressure gradients. Similarly, MCTS of a variety of tumor cell lines exhibit a 
concentric histo-morphology, with a necrotic core surrounded by a viable cell rim. The spherical symmetry 
is an important prerequisite for investigating the effect of environmental factors on cell proliferation and 
viability in a 3D environment on a quantitative basis [6]. 

2 Four-dimensional (4D) imaging of biological specimens (3D image reconstruction of the same living 
sample at different time points), is an application of confocal microscopy with fluorescence probes. 

3 Holographic coherence-domain imaging records full-frame depth resolved images throughout living 
multicellular tumor spheroids in vitro, without computed tomography, allowing real-time video fly-through 
under interactive control of the operator. 

4 Imaging bioluminescence technique allows the mapping of metabolite concentrations (e.g., ATP, glu- 
cose, and lactate) in cryosections of spheroid sections at a high spatial resolution [7]. 
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Figure 1: An example of multicellular tumor spheroid (MTS): a differential interference 
contrast image. 

parabolic reaction-diffusion PDeH 

= V • (A*i) - V • ( Vi $i) + Ai(*i, Ci) - nt($i, d) (1) 

(t?t = d/dt), where for phase i, $j is the volume fraction (Yli = 1)) Di is the random 
motility or diffusion, Aj($j, Cj) is the chemical and phase dependent production, and 
jUj($i,Ci) is the chemical and phase dependent degradation/death, and Vj is the cell 
velocity defined by the constitutive equation 

v< = -/xVp, (2) 

where /i is a positive constant describing the viscous-like properties of tumor cells and p 
is the spheroid internal pressure. 

In particular, the multi-phase equation (JT]) splits into two heat -like mass-conservation 
PDEs [3], 

d t $ c = S° -V -(v c $ c ), d t <S> F = S F - V-(v F $ F ), (3) 

where and are the tissue cell/matrix and fluid volume fractions, respectively, v c 
and are the cell/matrix and the fluid velocities (both defined by their constitutive 
equations of the form of ((2J)), S is the rate of production of solid phase tumor tissue and 
S F is the creation/degradation of the fluid phase. Conservation of matter in the tissue, 
$ C + $ F = 1, implies that V-(v c $° +v F <3? F ) = $ C + $ F . The assumption that the tumor 
may be described by two phases only implies that the new cell/matrix phase is formed 

5 Reaction-diffusion systems are PDE-models that describe how the concentration of one or more sub- 
stances distributed in space changes under the influence of two processes: local (bio) chemical reactions in 
which the substances are converted into each other, and diffusion which causes the substances to spread 
out in space. They have the form of semi-linear parabolic partial differential equations. 
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from the fluid phase and vice versa, so that S c + S F = 0. The detailed biochemistry of 
tumor growth can be coupled into the model above through the growth term S , with 
equations added for nutrient diffusion, see [3] and references therein. 

The multi-phase tumor growth model ([1]) has been derived from the classical trans- 
port/mass conservation equations for different chemical species inside the MTS [3j, 

d t d = Pi - V • Ni. (4) 

Here Cj are the concentrations of the chemical species, subindex a for oxygen, b for glucose, 
c for lactate ion, d for carbon dioxide, e for bicarbonate ion, / for chloride ion, and g for 
hydrogen ion concentration; Pj is the net rate of consumption/production of the chemical 
species both by tumor cells and due to the chemical reactions with other species; and Nj 
is the flux of each of the chemical species inside the tumor spheroid, given (in the simplest 
case of uncharged molecules of glucose, O2 and CO2) by Fick's law, 

Ni = -DiVCi, 

where D{ are (positive) constant diffusion coefficients. In case of charged molecules of 
ionic species, the flux Nj contains also the (negative) gradient of the volume fractions <3?j. 

In all above cases, tumor growth is, in terms of statistical mechanics, associated to 
entropy growth. The more uncertainty (measured as a number of tumor microstates) 
the tumor spheroid possesses, the larger is its entropy. Formally, we can apply Shannon 
formula to the probability distributions (p±, ■ ■ ■ ,p n ) of cancer cells within the body, 

S{pi,---,Pn) = -J2 Pil ° g2Pi - 

i 

In other words, if we do not control the tumor growth, naturally it is governed by the 
Second Law of Thermodynamics: 

d t S > 0, 

which is an expression of the universal law of increasing entropy, stating that the entropy 
(i.e., total number of tumor cells) of an isolated thermodynamic system (i.e., human body) 
which is not in equilibrium will tend to increase over time, approaching a maximum value 
at equilibrium (i.e., death threatening situation). 

On the other hand, the Ricci flow equation (or, the parabolic Einstein equation), in- 
troduced by R. Hamilton in 1982 [23], is the nonlinear heat -like evolution equation^ 

dtQij = ~2Rij, (5) 

6 The current hot topic in geometric topology is the Ricci flow, a Riemannian evolution machinery that 
recently allowed G. Perelman to prove the celebrated Poincare Conjecture, a century-old mathematics 
problem (and one of the seven Millennium Prize Problems of the Clay Mathematics Institute) - and win 
him the 2006 Fields Medal (which he declined in a public controversy) |14| . The Poincare Conjecture can 
roughly be put as a question: Is a closed 3D manifold M topologically a sphere if every closed curve in 
M can be shrunk continuously to a point? In other words, Poincare conjectured: A simply-connected 
compact 3D manifold is diffeomorphic to the 3D sphere 5* 3 (see e.g., [15]). 
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for a time-dependent Riemannian metric g = Qij{t) on a smooth reaQ n— manifold M with 
the Ricci curvature tensor Rij H This equation roughly says that we can deform any metric 
on a 2D surface or nD manifold by the negative of its curvature; after normalization (see 
Figure [2]), the final state of such deformation will be a metric with constant curvature. 
The factor of 2 in ([SJ is more or less arbitrary, but the negative sign is essential to insure 
a kind of complex volume exponential decays] since the Ricci flow equation ([5]) is a kind 
of nonlinear generalization of the standard linear heat equation 

dtu = An. (6) 

Like the heat equation (|6|), the Ricci flow equation ([5]) is well behaved in forward time 
and acts as a kind of smoothing operator (but is usually impossible to solve in backward 
time). If some parts of a solid object are hot and others are cold, then, under the heat 
equation, heat will flow from hot to cold, so that the object gradually attains a uniform 
temperature. To some extent the Ricci flow behaves similarly, so that the Ricci curvature 
'tries' to become more uniform |24j . thus depicting a monotonic entropy growth^^ dtS > 0, 



7 For the related Kahler Ricci flow on complex manifolds, see e.g., [161 117] 

8 This particular PDE ([5| was chosen by Hamilton for much the same reason that A. Einstein introduced 
the Ricci tensor into his gravitation field equation, 

Rij — 2^j^ = SnTij, 

where is the energy-momentum tensor. Einstein needed a symmetric 2-index tensor which arises 
naturally from the metric tensor gij and its first and second partial derivatives. The Ricci tensor Rij is 
essentially the only possibility. In gravitation theory and cosmology, the Ricci tensor has the volume- 
decreasing effect (i.e., convergence of neighboring geodesies, see [19|). 

9 This complex geometric process is globally similar to a generic exponential decay ODE: 

x = -\f(x), 

for a positive function f(x). We can get some insight into its solution from the simple exponential decay 
ODE, 

x = —\x with the solution x(t) = xoe~ Xt , 

(where x — x(t) is the observed quantity with its initial value xo and A is a positive decay constant), as 
well as the corresponding nth order rate equation (where n > 1 is an integer), 

11 

x = — \x n with the solution = + (n — 1) At. 

x n-i xq"- 1 

10 Note that two different kinds of entropy functional have been introduced into the theory of the Ricci 
flow, both motivated by concepts of entropy in thermodynamics, statistical mechanics and information 
theory. One is Hamilton's entropy, the other is Perelman's entropy. While in Hamilton's entropy, the 
scalar curvature R of the metric c?y is viewed as the leading quantity of the system and plays the role 
of a probability density, in Perelman's entropy the leading quantity describing the system is the metric 
gij itself. Hamilton established the monotonicity of his entropy along the volume-normalized Ricci flow 
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which is due to the positive definiteness of the metric gij > 0, and naturally implying the 
arrow of time [201 El EE] • 

In a suitable local coordinate system, the Ricci flow equation ([5]) has a nonlinear heat- 
type form, as follows. At any time t, we can choose local harmonic coordinates so that 
the coordinate functions are locally defined harmonic functions in the metric g(t). Then 
the Ricci flow takes the form (see e.g., |25| ) 

dtQij = &9ij + Qij(g, dg), (7) 

where A is the Laplace-Beltrami differential operator on functions with respect to the 
metric g and Q is a lower-order term quadratic in g and its first order partial deriva- 
tives. From the analysis of nonlinear heat PDEs, one obtains existence and uniqueness of 
forward-time solutions to the Ricci flow ([7]) on some time interval, starting at any smooth 
initial metric go = <?«(0). 

As a simple example of the Ricci flow equations ([5|)-([7|), consider a round spherical 
boundary S 2 of the MTS of radius r. The metric tensor on S 2 takes the form 

9ij = r Sij: 

where gij is the metric for a unit sphere, while the Ricci tensor 

Rij = (n - l)9ij 
is independent of r. The Ricci flow equation on S 2 reduces to 

r 2 = -2(n - 1), 

with solution 

r 2 (t) = r 2 (0) - 2(n - 

Thus the boundary sphere S 2 collapses to a point in finite time (see |24J). 

More generally, the following geometrization conjecture holds for an MTS 3-manifold 
M: Suppose that we start with a compact initial MTS-manifold Mq whose Ricci tensor 
Rij is everywhere positive definite. Then, as Mq shrinks to a point under the Ricci flow 
(|SJ), it becomes rounder and rounder. If we rescale the metric gij on Mq so that the volume 
of Mq remains constant, then Mq converges towards another compact MTS-manifold M\ 
of constant positive curvature (see [23]). 

In case of even more general MTS 3— manifolds (outside the class of positive Ricci 
curvature metrics), the situation is much more complicated, as various singularities may 
arise. One way in which singularities may arise during the Ricci flow is that a spherical 

on the 2-sphere S 2 [28] . Perelman established the monotonicity of his entropy along the Ricci flow in all 
dimensions |21j . 
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boundary S 2 = dM of an MTS 3-manifold M may collapse to a point in finite time. 
Such collapses can be eliminated by performing a kind of "geometric surgery" on the MTS 
manifold M, that is a sophisticated sequence of cutting and pasting without accumulation 
of time errorj^l (see [22]). After a finite number of such surgeries, each component either: 
(i) converges towards a 3-manifold of constant positive Ricci curvature which shrinks to a 
point in finite time, or possibly (ii) converges towards an S 2 x S 1 which shrinks to a circle 
S 1 in finite time, or (iii) admits a "thin-thick" decomposition of [30J. Therefore, one can 
choose the surgery parameters so that there is a well defined Ricci-flow-with surgery, that 
exists for all time |22j . 

In this paper we use the evolving geometric machinery of the volume-decaying and 
entropy-growing Ricci flow g(t), given by equations ([5])-([7|), for modelling general 3D 
avascular MTS decay, corresponding to parabolic multi-phase reaction-diffusion PDEs 



2 Ricci flow and multi-phase avascular MTS decay control 
2.1 Geometrization Conjecture 

Recall that geometry and topology of smooth surfaces are related by the Gauss-Bonnet 
formula for a closed surface S (see, e.g., [2| [T7]) 



where dA is the area element of a metric g on E, K is the Gaussian curvature, x(E) is 
the Euler characteristic of E and gen(E) is its genus, or number of handles, of E. Every 
closed surface E admits a metric of constant Gaussian curvature K = +1, 0, or —1 and 
so is uniformized by elliptic, Euclidean, or hyperbolic geometry, which respectively have 
gen(5" 2 ) = (sphere), gen(T 2 ) = 1 (torus) and gen(E) > 1 (torus with several holes). The 
integral ([8]) is a topological invariant of the surface E, always equal to 2 for all topological 
spheres S 2 (that is, for all closed surfaces without holes that can be continuously deformed 
from the geometrical sphere) and always equal to for the topological torus T 2 (i.e., for 
all closed surfaces with one hole or handle). 

The general topological framework for the Ricci flow ([5]) is Thurston's Geometrization 
Conjecture [30], which states that the interior of any compact 3-manifold can be split in 

11 Hamilton's idea was to perform surgery to cut off the singularities and continue his flow after the 
surgery. If the flow develops singularities again, one repeats the process of performing surgery and contin- 
uing the flow. If one can prove there are only a finite number of surgeries in any finite time interval, and 
if the long-time behavior of solutions of the Ricci flow © with surgery is well understood, then one would 
be able to recognize the topological structure of the initial manifold. Thus Hamilton's program, when 
carried out successfully, would lead to a proof of the Poincare conjecture and Thurston's geometrization 
conjecture [15] . 



(HMD. 




(8) 
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an essentially unique way by disjoint embedded 2D spheres S 2 and tori T 2 into pieces 
and each piece admits one of 8 geometric structures (including (i) the 3D sphere S s with 
constant curvature +1; (ii) the 3D Euclidean space M 3 with constant curvature and (hi) 
the 3D hyperbolic space H 3 with constant curvature — 1)H The geometrization conjecture 
(which has the Poincare Conjecture as a special case) would give us a link between the 
geometry and topology of MTS 3-manifolds, analogous in spirit to the case of 2D surfaces. 

In higher dimensions, the Gaussian curvature K corresponds to the Riemann curvature 
tensor 9\m on a smooth n— manifold M, which is in local coordinates on M denoted by 
its (4, 0)— components R%jkh or its (3, 1)— components (see Appendix, as well as e.g., 
[21 [T7]). The trace (or, contraction) of 9im, using the inverse metric tensor g %1 = (gij) , is 
the Ricci tensor IKc, the 3D curvature tensor, which is in a local coordinate system {x l }'^ =1 
defined in an open set U C M, given by 

Rij = tr(«Rm) = g kl R ijk i 

(using Einstein's summation convention), while the scalar curvature is now given by the 
second contraction of 9\m as 

R = ix(p\c) = g ij Rij. 

In general, the Ricci flow gij(t) is a one-parameter family of Riemannian metrics on 
a compact n— manifold M governed by the equation ([5]), which has a unique solution for 
a short time for an arbitrary smooth metric g^ on M [23]. If 9\c > at any local point 
x = {x 1 } on M, then the Ricci flow ([5]) contracts the metric gij(t) near x, to the future, 
while if yic < 0, then the flow ([5]) expands gij(t) near x. The solution metric gijit) of the 
Ricci flow equation ([5]) shrinks in positive Ricci curvature direction while it expands in 
the negative Ricci curvature direction, because of the minus sign in the front of the Ricci 
tensor R^. In particular, in 2D, on a sphere S 2 , any metric of positive Gaussian curvature 
will shrink to a point in finite time. At a general point, there will be directions of positive 
and negative Ricci curvature along which the metric will locally contract or expand (see 
|25j). In 3D, if a simply-connected compact 3-manifold M has a Riemannian metric gij 
with positive Ricci curvature then it is diffeomorphic to the 3-sphere 5 3 |23j . 

All three Riemannian curvatures (R, ^Hc and 9tm), as well as the associated volume 
forms, evolve during the Ricci flow ([5]). 



12 Another five allowed geometric structures are represented by the following examples: (iv) the product 

S 2 x S 1 ; (v) the product M 2 x S 1 of hyperbolic plane and circle; (vi) a left invariant Riemannian metric on 

the special linear group SL(2, R); (vii) a left invariant Riemannian metric on the solvable Poincare-Lorentz 

group £7(1, 1), which consists of rigid motions of a (1 + 1)— dimensional space-time provided with the fiat 

metric dt 2 — dx 2 ; (viii) a left invariant metric on the nilpotent Heisenberg group, consisting of 3 x 3 matrices 

of the form 
1 * * 

1 * .In each case, the universal covering of the indicated manifold provides a canonical model 
1 J 
for the corresponding geometry [24] . 
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2.2 MTS evolution under the Ricci flow 



The Ricci flow evolution equation for the metric tensor gij implies the evolution equa- 
tion for the Riemann curvature tensor 9\m, 

dtVlm = Amm + Q n , (9) 

where Q n is a certain quadratic expression of the Riemann curvatures. From the general 
n— curvature expression ([9]) we have two special cases important for MTS-evolutionf^l 



• The 3D evolution equation for the Ricci curvature tensor 9\c on an MTS 3-manifold 

M, 

d t ftc = A<Kc + Q 3 , (10) 
where Q3 is a certain quadratic expression of the Ricci curvatures; and 

• The 2D evolution equation for the scalar surface curvature R, 

d t R = AR + 2\9\c\ 2 , (11) 

which holds both on an MTS 3-manifold M and on its 2D boundary surface dM. 
Therefore, by the maximum principle, the minimum of R is non-decreasing along 
the flow g(t), both on M and on dM (see [2"T]). 

Let us now see in detail how various MTS-related geometric quantities evolve given 
the short-time solution of the Ricci flow equation ([5]) on an MTS 3-manifold M. Let us 
first calculate the variation formulas for the Christoffel symbols and curvature tensors 
on M and then the corresponding evolution equations (see [231 ED E2])- If g(s) is a 
one-parameter family of metrics on M with 

13 By expanding the maximum principle for tensors, Hamilton proved that Ricci flow g(t) given by (JS| 
preserves the positivity of the Ricci tensor JHc in 3D (as well as of the Riemann curvature tensor SHm in all 
dimensions); moreover, the eigenvalues of the Ricci tensor in 3D (and of the curvature operator 5Hm in 4D) 
are getting pinched point- wisely as the curvature is getting large 23, 26 . This observation allowed him 
to prove the convergence results: the evolving metrics (on a compact manifold) of positive Ricci curvature 
in 3D (or positive Riemann curvature in 4D) converge, modulo scaling, to metrics of constant positive 
curvature. 

However, without assumptions on curvature, the long time behavior of the metric evolving by Ricci flow 
may be more complicated [3T|. In particular, as t approaches some finite time T, the curvatures may 
become arbitrarily large in some region while staying bounded in its complement. On the other hand, 
Hamilton 27 discovered a remarkable property of solutions with nonnegative curvature tensor iHm in 
arbitrary dimension, called the differential Harnack inequality, which allows, in particular, to compare the 
curvatures of the solution of © at different points and different times. 
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then the variation of the Christoffel symbols T k j on M is given by 

9,1* = -g kl (ViVji + VjVu ~ ViVij) , (12) 

from which follows the evolution of the Christoffel symbols I* under the Ricci flow g(t) 
on M given by (|5|), 

dtT^ = -g kl (VilJ,-, + VjRu - ViRti) . 
From (|12p we calculate the variation of the Ricci tensor Rij on M as 

d a Rij = V m (d s T%) - V< (5 s r™ ) , (13) 

and the variation of scalar curvature R on M by 

d s R = -AF + div(divu) - (v,mc), (14) 

where V = g^Vy = tr(-u) is the trace of v = (%•)• 

If an MTS 3— manifold M is oriented, then the volume 3— form on M is given, in a 
positively oriented local coordinate system {x % } E U C M, b}0 

dfi = J det(gij) dx 1 A cte 2 A dx 3 . (15) 

If d s #ij = %, then 

The evolution of the volume form d/i under the Ricci flow g(t) on M is given by the 
exponential decay/growth relation with the scalar curvature R as the rate constant, 

dfdfj, = —Rdfi, (16) 

which gives an exponential decay for R > (elliptic geometry) and exponential growth for 
R < (hyperbolic geometry). The elementary volume evolution (|16p implies the integral 
form of the exponential relation for the total MTS-volume 

vol (3) = / d ^ 

J M 



in the form 



d t volGK*)) = - / Rdfi, 

<M 



14 Extension to higher-dimensional Riemannian manifolds is obvious |17] : also, for related volume forms 
on symplectic manifolds, see [18] 
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which again gives an exponential decay for elliptic R > and exponential growth for 
hyperbolic R < 0. 

This is a crucial point for the tumor decay control: we need to keep the elliptic geometry 
of the MTS - by all possible means. And naturally - it will be so, because it started as a 
spherical shape with R > 0. We just need to keep the MTS in this shape and prevent any 
hyperbolic distortions of R < 0. The, it will naturally have an exponential decay. 

On the other hand, if we are not able to keep the positivity of the scalar curvature of 
the MTS, and thus prevent its expanding to in infinity ('deadly' hyperbolic case), we can 
also consider the normalized Ricci flow of the MTS on its 3-manifold M (see Figure [2] as 
well as ref. [3T]): 

2 

dtfjij = -2Rij + -mj, (17) 

where 

f = vol(^) -1 / Rdfi 
Jm 

is the average scalar curvature on M. We then have the MTS volume conservation law: 

d t vo\(g(t)) = 0. 

♦ t • 

Figure 2: An example of Ricci flow normalization: unnormalized flow (up) and normalized 
flow (down). 

To study the long-time existence of the normalized Ricci flow (|17p on an MTS 3— mani- 
fold M, it is important to know what kind of curvature conditions are preserved under the 
equation. In general, the Ricci flow g(t) on M tends to preserve some kind of positivity 
of curvatures. For example, positive scalar curvature R is preserved both on M and 
on its boundary dM. This follows from applying the maximum principle to the evolution 
equation (jlip for scalar curvature R both on M and on dM. Also, positive Ricci curvature 
is preserved under the Ricci flow on M. (This is a special feature of 3D and is related to 
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the fact that the Riemann curvature tensor may be recovered algebraically from the Ricci 
tensor and the metric in 3D |31j.) 

In particular, we have the following result (see [28]) for MTS-surfaces dM: Let dM 
be a closed MTS-surface. Then for any initial 2D metric go on dM, the solution to 
the normalized Ricci flow (|17p on dM exists for all time. Moreover, (i) If the Euler 
characteristic of dM is non-positive, then the solution metric g(t) on dM converges to 
a constant curvature metric as t — ► oo; and (ii) If the scalar curvature R of the initial 
metric go is positive, then the solution metric g(t) on dM converges to a positive constant 
curvature metric as t — > oo. (For surfaces with non-positive Euler characteristic, the proof 
is based primarily on maximum principle estimates for the scalar curvature.) 

In other words, the normalized Ricci flow of the MTS will make it completely round 
with a geometrical sphere shell - ideal for surgical removal. This is our second option for 
the MTS control. If we cannot force it to exponential decay, then we must try to normalize 
into a round spherical shell - which is suitable for surgical removal. 

The negative flow of the total MTS-volume vo\(g(t)) is the Einstein-Hilbert functional, 
given by (see 03J ED [25] ) 

E(g) = [ Rd» = -d t vol(g(t)). 

If we put d s gij = Vij , we have 

d s E{g) = J (- AV + div(div v) - (v, 9\c) + \RV^J dfi 

= j M (v, ^Rgij - Rij^ dfi, 

so the critical points of E(g) satisfy Einstein's equation 

1 

2^9ij ~ Rij — 0. 

The gradient flow of E(g) on M, given by 

d t9ij = 2 (VB( ff )).. = Rg i3 - 2R ij , 

is almost the Ricci flow (JSJ. Thus, Einstein metrics are the fixed points of the Ricci flow 
g(t) on M0 

15 In 3D manifolds, Einstein metrics are metrics with constant curvature. However, along the way, the 
deformation will encounter singularities. The major question, resolved by Perelman, was how to find a 
way to describe all possible singularities. 
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Let A denote the Laplacian acting on functions on an MTS 3— mani-fold M, which is 
in local coordinates {x 1 } € U C M given by 



A = gVViVj = g ij (dij - J%dk) . 



For any smooth function / on M we have (23^ [32] 



A|V/| 2 



ViA/ + RijVjf, and 

2|V i V i /| 2 + 2R ij V l fV J f + 2V i /V l A/. 



From this it follows that if we have 



VV/ = 



<Kc > 



and 



A/^0 



), |V/| = 1, 
5Hc(V/,V/) = 0. 



then 



Using A, we can write the linear heat equation on M as 



dtu = Au 



where u is the MTS-temperature. In particular, the Laplacian acting on functions with 
respect to g(t) will be denoted by A g r t \. If (M, g(t)) is a solution to the Ricci flow equation 
([5]), then we have 



Now, the evolution equation (fTT|) for the scalar curvature R under the Ricci flow © 
follows from ()14|) . Using equation (|27p from Appendix, we have: 



showing again that the scalar curvature R satisfies a heat-type equation with a quadratic 
nonlinearity both on an MTS 3-manifold M and on its 2D boundary surface dM. 

Next we will find the exact form of the evolution equation (|1U|) for the Ricci tensor £Hc 
under the Ricci flow g(t) given by © on an MTS 3-manifold M. (Note that in higher 
dimensions, the appropriate formula would involve the whole Riemann curvature tensor 
9\m.) In general, given a variation d s gij = Vij, from (JT3J) we get 



where A^ denotes the so-called Lichnerowicz Laplacian (which depends on 9\m) (see 




— 2RijViVj. 





[231 [32]). Since 



ViVjR - V f (div(JRc))j - Vj(div(9tc))i = 0, 
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by (|27|) (after some algebra) we get that under the Ricci flow © the evolution equation 
for the Ricci tensor 9\c on M is 

dtRij = Ai?jj + 3RRij — QRimRjm + (2|£Kc| 2 — i? 2 ) gij. 

So, just as in case of the evolution ([! ip of the scalar curvature dtR (both on M and on 
its boundary dM), we get a heat-type evolution equation with a quadratic nonlinearity 
for dtRij, which means that positive Ricci curvature (9tc > 0) of elliptic MTS-geometry 
is preserved under the Ricci flow g{t) on M. 

More generally, we have the following result for MTS 3-manifolds (see [23]): Let (M, go) 
be a compact Riemannian MTS 3— manifold with positive Ricci curvature £Hc. Then there 
exists a unique solution to the normalized Ricci flow g(t) on M with g(0) = go for all time 
and the metrics g(t) converge exponentially fast to a constant positive sectional curvature 
metric g^ on M. In particular, M is diffeomorphic to a 3D sphere 5 3 . (As a consequence, 
such an MST 3— manifold M is necessarily diffeomorphic to a quotient of the 3— sphere by 
a finite group of isometries. It follows that given any homotopy 3— sphere, if one can show 
that it admits a metric with positive Ricci curvature, then the Poincare Conjecture would 
follow [31].) In particular, compact and closed 3— manifolds which admit a non-singular 
solution can also be decomposed into geometric pieces [29] . 



2.3 Ricci breathers and solitons 

Recall that breathers are solitonic structures given by localized periodic solutions of some 
nonlinear soliton PDEs, including the exactly solvable sine-Gordon equation^! and the 
focusing nonlinear Schrodinger equation 1^1 



Hi 



An exact solution u = u(x,t) of the (1+1)D sine-Gordon equation 

d 2 u _ d 2 u 
dt 2 dx 2 

is J37] 

y/1 — LO 2 COs(drf) 



u = 4 arctan 



yio cosh(\/l — LO 2 x) , 

which, for u < 1, is periodic in time t and decays exponentially when moving away from x — 0. 
17 The focusing nonlinear Schrodinger equation is the dispersive complex-valued (1+1)D PDE [38] , 

. du d 2 u . l2 _ , . / — 7 , xx 

* ^7 + + \ u \ u = °> (* = V -1, u = u[x,t)) 
at ox A 

with a breather solution of the form: 

( 2b 2 cosh(6>) + 2 ib V2=¥ sinh(g) \ 2 2 /— — 

u = — ■= — , — — 1 a expfj a t) with & = a b\J2 — b i t, 

V 2 cosh(6») - V2 V2 - & 2 cos(abs) / 

which gives breathers periodic in space x and approaching the uniform value a when moving away from 
the focus time t = 0. 
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A metric gij{t) evolving by the Ricci flow g(t) given by ([5]) on an MTS 3-manifold M 
is called a Ricci breather, if for some t\ < t 2 and a > the metrics agij(t\) and gij(t 2 ) 
differ only by a diffeomorphism; the cases a = l,a < l,a > 1 correspond to steady, 
shrinking and expanding breathers, respectively. Trivial breathers on M, for which the 
metrics <ft_j(ti) and gij(t 2 ) differ only by diffeomorphism and scaling for each pair of t\ 
and t 2 , are called Ricci solitons. Thus, if one considers Ricci flow as a dynamical system 
on the space of Riemannian metrics modulo diffeomorphism and scaling, then breathers 
and solitons correspond to periodic orbits and fixed points respectively. At each time the 
Ricci soliton metric satisfies on M an equation of the form [21] 

Rij + cgij + Vibj + Vjbi = 0, 

where c is a number and 6j is a one-form; in particular, when bi = ^V^a for some function 
a on M, we get a gradient Ricci soliton. An important example of a gradient shrinking 
soliton is the Gaussian soliton, for which the metric gij is just the Euclidean metric on 
R 3 , c = 1 and a = -\x\ 2 /2. 

2.4 Heat equation and Ricci entropy 

Given a C 2 function u : M — > R on a Riemannian MTS 3— manifold M, its Laplacian is 
defined in local coordinates {x 1 } G U C M to be 

An = tig (V 2 n) = g ij ViVjU, 

where Vj is its associated covariant derivative (Levi-Civita connection, see Appendix). 
We say that a C 2 function u : M x [0, T) — > IR, where T £ (0, oo], is a solution to the heat 
equation on M if 

dtu = Au. (18) 

One of the most important properties satisfied by the heat equation is the maximum 
principle, which says that for any smooth solution to the heat equation, whatever point- 
wise bounds hold at t = also hold for t > |31j. More precisely, we can state: Let 
u : M x [0, T) — > R be a C 2 solution to the heat equation (|18p on a complete Riemannian 
MTS 3-manifold M. H d < u {x, 0) < C 2 for all x £ M, for some constants C u C 2 G R, 
then Ci < it (x, t) < C2 for all x G M and t G [0, T). This property exhibits the smoothing 
behavior of the heat equation (|18p on M. 

Now, consider Perelman's entropy functional [21] on an MTS 3-manifold M 

?= f (R+\Vf\ 2 )e- f dfx (19) 
Jm 

for a Riemannian metric g™ and a (temperature-like) scalar function / on a closed 3— 
manifold M, where dfi is the volume 3-form (I15p . During the Ricci flow ([5]), J 7 evolves on 
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M as 



d t T = 2 J {Rij + V i V J /| 2 e- / d A t. (20) 
Now, define \{g%j) = inf !F(gij, /), where infimum is taken over all smooth /, satisfying 

e -/ d/u = l. (21) 



M 



\(gij) is the lowest eigenvalue of the operator — 4A+i?. Then the entropy evolution formula 
([20]) implies that X(gij(t)) is nondecreasing in t, and moreover, if X(h) = Xfo), then for 
t £ [ii, £2] we have + VjVj/ = for / which minimizes f onM |21j . Thus a steady 
breather on M is necessarily a steady soliton. 

If we define the conjugate heat operator on M as 

□* = -d/dt -A + R 

then we have the conjugate heat equatio 



D*u = 0. (23) 



18 In [21] Perelman stated a differential Li-Yau-Hamilton (LYH) type inequality [33] for the fundamental 
solution it = u(x,t) of the conjugate heat equation (|23[) on a closed n— manifold M evolving by the Ricci 



flow (HI). Let p G M and 



u = (47rr) 2 e f 



be the fundamental solution of the conjugate heat equation in M x (0, T), 

□*u = 0, or dtu + Ait = Ru, 

where r = T — t and R = R(-,t) is the scalar curvature of M with respect to the metric g(t) with 
lim t u — 5 P (in the distribution sense), where 5 P is the delta-mass at p. Let 

v = [r(2A/ - I V/| 2 +#)+/- n]w, 

where r = T — i. Then we have a differential LYH-type inequality 

v(x,t) < in M x (0,T). (22) 

This result was used by Perelman to give a proof of the pseudolocality theorem [21] which roughly said 
that almost Euclidean regions of large curvature in closed manifold with metric evolving by Ricci flow g(t) 
given by (J5j remain localized. 

In particular, let (M,g(t)), < t < T, dM 7^ (f>, be a compact 3— manifold (like MTS) with metric g(t) 
evolving by the Ricci flow g(t) given by (0 such that the second fundamental form of the surface dM with 
respect to the unit outward normal d/du of dM is uniformly bounded below on dM x [0,T]. A global 
Li-Yau gradient estimate [34] for the solution of the generalized conjugate heat equation was proved in 
|33| (using a a variation of the method of P. Li and S.T. Yau, [31]) on such a manifold with Neumann 
boundary condition. 
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The entropy functional (fT9j) is nondecreasing under the following coupled Ricci-heat 
flow on M [35] 

9t9ij = — 2i?jj, 

d t u = -Au -J L + -w, 24 

li 2 

where the modified conjugate heat equation (|24|) ensures 



K 2 d/i, = 1 

M 

f 

to be preserved by the Ricci flow g(t) on M. If we define u = e~ 2 , then (|24^) is equivalent 
to /—evolution equation on M, 

9 t / = -A/+|V/| 2 -i?, 

which instead preserves (|2ip. 

2.5 Thermodynamic analogy 

Perelman's functional is analogous to negative entropy [21]. Recall that thermodynamic 
partition function for a generic canonical ensemble at temperature (3~ l is given by 

Z = j e- 0E duj( E ), (25) 

where u(E) is a 'density measure', which does not depend on (3. From it, the average 
energy is given by 

(E) = -d lnZ, 

the entropy is 

S = /3(E) + In Z, 

and the fluctuation is 

<x=((£-(£)) 2 }=^ 2 lnZ. 

If we now fix a closed MTS 3— manifold M with a probability measure m and a metric 
gij{r) that depends on the temperature r, then according to equation 

d T g lj = 2(R ij +V i V J f), 

the partition function (I25p is given by 

lnZ= /(-/ + |)dm. (26) 
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Prom ([26]) we get (see [21]) 



(E) 



—T 



2 / (i?+|V/| 2 -f )dm, 
Jm It 



S 



[ (r(R+\Vf\ 2 ) + f-n)dm, 



Jm 




Jm 



[ \Rij + VjV jf - ^-gij\ 2 dm 



where 



dm = udV, 



u = (Attt) 



From the above formulas, we see that the MTS-fluctuation a is nonnegative; it vanishes 
only on a gradient shrinking soliton. (E) is nonnegative as well, whenever the flow exists 
for all sufficiently small r > 0. Furthermore, if the MTS-heat function u: (a) tends to a 
S— function as r — > 0, or (b) is a limit of a sequence of partial heat functions Uj, such that 
each Ui tends to a 5— function as r — > n > 0, and r» — > 0, then the MTS-entropy S is 
also nonnegative. In case (a), all the quantities (E) ,S,a tend to zero as r — > 0, while in 
case (b), which may be interesting if gij(r) becomes singular at r = 0, the MTS-entropy 
S may tend to a positive limit. 

3 Monoclonal antibodies for MTS— decay control 

To keep the MTS within the elliptic geometry with the positive scalar curvature, R > 0, 
which would enable the exponential decay of its volume, we need the help from the local 
immune system. 

Recall that monoclonal antibodies (mAb) are monospecific antibodies that are identical 
because they are produced by one type of immune cell that are all clones of a single parent 
cell. Given (almost) any substance, it is possible to create monoclonal antibodies that 
specifically bind to that substance; they can then serve to detect or purify that substance 



The idea of a 'magic bullet' was first proposed by Paul Ehrlich (Nobel Prize in Phys- 
iology or Medicine in 1908), who a century a go postulated that if a compound could be 
made that selectively targeted a disease-causing organism, then a toxin for that organism 
could be delivered along with the agent of selectivity. 

The invention of monoclonal antibodies is generally accredited to Georges Khler, Csar 
Milstein, and Niels Kaj Jerne in 1975 |40j . who shared the Nobel Prize in Physiology or 
Medicine in 1984 for the discovery. The key idea was to use a line of myeloma cells that 
had lost their ability to secrete antibodies, come up with a technique to fuse these cells 
with healthy antibody producing B-ceils, and be able to select for the successfully fused 
cells. 



[39j. 
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Human monoclonal antibodies are produced using transgenic mice or phage display 
libraries. Human monoclonal antibodies are produced by transferring human immunoglob- 
ulin genes into the murine genome, after which the transgenic mouse is vaccinated against 
the desired antigen, leading to the production of monoclonal antibodies. Phage display 
libraries allow the transformation of murine antibodies in vitro into fully human antibod- 
ies. 

Antibody-directed enzyme prodrug therapy (ADEPT) involves the application of can- 
cer associated monoclonal antibodies which are linked to a drug-activating enzyme. Sub- 
sequent systemic administration of a non-toxic agent results in its conversion to a toxic 
drug, and resulting in a cytotoxic effect which can be targeted at malignant cells. The 
clinical success of ADEPT treatments has been limited to date |41j . However it holds 
great promise, and recent reports suggest that it will have a role in future oncological 
treatment |42j . 

4 Appendix: Riemann and Ricci curvatures on a smooth 
n— manifold 

Recall that proper differentiation of vector and tensor fields on a smooth Riemannian 
n— manifold is performed using the Levi-Civita covariant derivative (see, e.g., [21 [T7]). 
Formally, let M be a Riemannian n— manifold with the tangent bundle TM and a local 
coordinate system {x l }^ =1 defined in an open set U C M. The covariant derivative 
operator, V x ■ C°°(TAf) — > C°°(TM), is the unique linear map such that for any vector 
fields X, Y, Z, constant c, and function / the following properties are valid: 

Vx+cY = Vx+cVy, 

V x (Y + fZ) = V x Y + (Xf)Z + fV x Z, 
VxY-VyX = [X,Y], 

where [X, Y] is the Lie bracket of X and Y (see, e.g., [18J). In local coordinates, the metric 
g is defined for any orthonormal basis (pi = d x %) in U C M by 

9ij = g(di,dj) = Sij, d k gij = 0. 

Then the affine Levi-Civita connection is defined on M by 

Vdpj = Tijdk, where 

are the (second-order) Christoffel symbols. 
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Now, using the covariant derivative operator Vx we can define the Riemann curvature 
(3, l)-tensor <Km by (see, e.g., 0Q7]) 

mm(X, Y)Z = V x VyZ - V Y V X Z - V [x ,y] z - 

9\m measures the curvature of the manifold by expressing how noncommutative covariant 
differentiation is. The (3, 1)— components R\- k of £Km are defined in U C M by 

9\m(di,dj)dk = R^di, which expands as [35], 

r>l f\ -pi f\ -pi _i_ -pm-pl pm-pi 

^Hjk ~ u i L jk u j L ik ' L jk L im 1 ik 1 jm m 

Also, the Riemann (4,0)— tensor Rijki = QimR^jk ls defined as the g— based inner product 
on M, 

Rijki = (^tn (pi, dj) d k , di) . 
The first and second Bianchi identities for the Riemann (4,0)— tensor Rijki hold, 

Rijki + Rjkil + Rkijl = 0, 
^iRjklm + V ' jRkilm + V 'kRijlm = 0, 

while the twice contracted second Bianchi identity reads 

2V j R lj = ViR. (27) 
The (0,2) Ricci tensor 9\c is the trace of the Riemann (3, 1)— tensor D^m, 

<Rc(y, Z) + tr(X -» fHm(X,r)Z), 

so that 

9tc(x,y) = 5 ($nm(a 4 ,x)a l ,y), 

Its components Rjk = (dj,dk) are given in t7 C M by the contraction (see e.g., [33]) 
i^jfc = or, in terms of Christoffel symbols, 

p pi pi a pi I pi ptn _ pi pra 

•ryfc — u i L jk u k L ji ~r L mi L jk L mk L ji- 

Being a symmetric second-order tensor, d\c has (p^-) independent components on an 
n— manifold M. In particular, on a 3-manifold, it has 6 components, and on a 2D surface 
it has only the following 3 components: 

-Rn = g 22 R2U2, R12 = g 12 R2m, R22 = g n Ri22i, 
which are all proportional to the corresponding coordinates of the metric tensor, 

R11 _ R12 _ R22 _ -R1212 ^g) 
5ii 512 522 det(g)' 
Finally, the scalar curvature R is the trace of the Ricci tensor JHc, given in U C M by 

R = g ij R i3 . 



21 



References 



[1] R.A. Gatenby, P.K. Maini, Mathematical oncology: Cancer summed up. Nature, 421, 
321, (2003). 

[2] V. Ivancevic, T. Ivancevic, Geometrical Dynamics of Complex Systems. Springer, 
Dordrecht, (2006). 

[3] T. Roose, S.J. Chapman, P.K. Maini, Mathematical Models of Avascular Tumor 
Growth. SIAM Rev. 49(2), 179-208, (2007). 

[4] L.A. Kunz-Schughart, Multicellular tumor spheroids: intermediates between mono- 
layer culture and in vivo tumor, Cell Biol. Int. 23(3), 157-61, (1999). 

[5] R.M. Sutherland, Cell and environment interactions in tumor microregions: The 
multicell spheroid model, Science, 240, 177-184, (1988). 

[6] S. Walenta, J. Doetsch, W. Mueller-Klieser, L.A. Kunz-Schughart, Metabolic Imag- 
ing in Multicellular Spheroids of Oncogene-transfected Fibroblasts, J. Histochem. 
Cythochem. 48, 509-522, (2000). 

[7] S. Walenta, J. Doetsch, W. Mueller-Klieser, ATP concentrations in multicellular 
spheroids assessed by single photon imaging and quantitative bioluminescence. Eur. 
J. Cell. Biol. 52, 389-393, (1990). 

[8] L. Preziosi, Cancer Modeling and Simulation. CRC Press, (2003). 

[9] G. Hu, Dongqing Li, Three-dimensional modeling of transport of nutrients for multi- 
cellular tumor spheroid culture in a microchanne, Biomed. Microdev. 9(3), 315-323, 
(2007). 

[10] S. Khoshyomn, P.L. Penar, W.J. McBride, D.J. Taatjes, Four-dimensional analysis of 
human brain tumor spheroid invasion into fetal rat brain aggregates using confocal 
scanning laser microscopy, J. Neuro-Oncol., 38(1), 1-10, (1998). 

[11] K. Jeong, Fourier-domain holographic optical coherence imaging of tumor spheroids 
and mouse eye, Appl. Opt. 44, 1798-1806, (2005). 

[12] T. Ivancevic, L. Jain, J. Pattison, A. Hariz, Nonlinear Dynamics and Chaos Methods 
in Neurodynamics and Complex Data Analysis. Nonl. Dyn. (Springer) (to appear) 

[13] V. Ivancevic, T. Ivancevic, High-Dimensional Chaotic and Attractor Systems. 
Springer, Berlin, (2006). 



22 



[14] D. Mackenzie, Perelman Declines Math's Top Prize; Three Others Honored in Madrid, 
Science 313, 1027, (2006). 

[15] S.T. Yau, Structure of Three-Manifolds - Poincare and geometrization conjectures, 
talk given at the Morningside Center of Mathematics on June 20, (2006). 

[16] V. Ivancevic, T. Ivancevic, Complex Dynamics: Advanced System Dynamics in Com- 
plex Variables. Springer, Dordrecht, (2007) 

[17] Ivancevic, V., Ivancevic, T.: Applied Differfential Geometry: A Modern Introduction. 
World Scientific, Singapore, (2007) 

[18] V. Ivancevic, Symplectic Rotational Geometry in Human Biomechanics, SIAM Rev. 
46(3), 455-474, (2004) 

[19] S. Hawking, R. Penrose, The Nature of Space and Time, Princeton Univ. Press, 1996. 

[20] R. Penrose, Singularities and Time-Asymmetry, in General Relativity: An Einstein 
Centenary Survey, (ed. S. Hawking, W. Israel), 581-638, Cambridge Univ. Press, 
(1979). 

[21] G. Perelman, The entropy formula for the Ricci flow and its geometric applications, 
arXiv:math.DG/0211159 

[22] G. Perelman, Ricci flow with surgery on three-manifolds, arXiv:math.DG/0303109 

[23] R.S. Hamilton, Three-manifolds with positive Ricci curvature, J. Diff. Geom. 17, 
255-306, (1982) 

[24] J. Milnor, Towards the Poincare Conjecture and the Classification of 3-Manifolds, 
Not. Am. Math. Soc. 50(10), 1226-1233, (2003) 

[25] M.T. Anderson, Geometrization of 3-manifolds via the Ricci flow, Not. Am. Math. 
Soc. 51(2), 184-193, (2004) 

[26] R.S. Hamilton, Four-manifolds with positive curvature operator, J. Dif. Geom. 24 
(1986) 153-179. 

[27] R.S. Hamilton, The Harnack estimate for the Ricci flow, J. Dif. Geom. 37 (1993) 
225-243. 

[28] R.S. Hamilton, The Ricci flow on surfaces, Cont. Math. 71 (1988), 237-261. 

[29] R.S. Hamilton, Non-singular solutions of the Ricci flow on three-manifolds, Comm. 
Anal. Geom., 7(4), 695-729, 1999. 



23 



[30] W. Thurston, Three-dimensional manifolds, Kleinian groups and hyperbolic geome- 
try, Bull. Amer. Math. Soc. 6 (1982) 357-381. 

[31] H.D., Cao, B. Chow, Recent developments on the Ricci flow, Bull. Amer. Math. Soc. 
36 59-74, 1999. 

[32] Chow, B., Knopf, D. The Ricci flow: An introduction, Mathematical Surveys and 
Monographs, AMS, Providence, RI, 2004. 

[33] S.Y. Hsu, Some results for the Perelman LYH-type inequality, 
arXiv:math.DG/0801.3506 

[34] P. Li, S.T. Yau, On the parabolic kernel of the Schrodinger operator. Acta Math. 
156, 1986, 153-201. 

[35] R. Ye, The log entropy functional along the Ricci flow. arXiv:math.DG/0708.2008v3 

[36] J. Li, First variation of the Log Entropy functional along the Ricci flow, 
arXiv:math.DG /0712.0832 

[37] M.J. Ablowitz, D.J. Kaup, A.C. Newell, H. Segur, Method for solving the sine-Gordon 
equation. Phys. Rev. Let. 30, 1262-1264, (1973). 

[38] N.N. Akhmediev, V.M. Eleonskii, N.E. Kulagin, First-order exact solutions of the 
nonlinear Schrodinger equation. Th. Math. Physics 72, 809-818, (1987). 

[39] Wikipedia, the free encyclopedia, 2008. 

[40] Kohler G, Milstein C, Continuous cultures of fused cells secreting antibody of pre- 
defined specificity. Nature, 256, 495, (1975). 

[41] Francis RJ, Sharma SK, Springer C, et al. A phase I trial of antibody directed enzyme 
prodrug therapy (ADEPT) in patients with advanced colorectal carcinoma or other 
CEA producing tumours. Br J Cancer 2002; 87: 600-607. 

[42] Carter P, Improving the efficacy of antibody-based cancer therapies. Nat. Rev. Can- 
cer, 1, 118-129, (2001). 

[43] C. Misner, K. Thorne, J. A. Wheeler, Gravitation, W.H. Freeman and Company, 
(1973). 



24 



